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Abstract. In anisotropic or layered superconductors thermal fluctuations as well as impurities 
induce a van der Waals (vdW) attraction between flux lines, as has recently been shown by Blatter 
and Geshkenbein in the thermal case [Phys. Rev. Lett. 77, 4^958 (1996)] and by Mukherji and 
Nattermann in the disorder dominated case [Phys. Rev. Lett. 79, 139 (1997)]. This attraction 
together with the entropic or disorder induced repulsion has interesting consequences for the low 
field phase diagram. We present two derivations of the vdW attraction, one of which is based 
on an intuitive picture, the other one following from a systematic expansion of the free energy 
of two interacting flux lines. Both the thermal and the disorder dominated case are considered. 
In the thermal case in the absence of disorder, we use scaling arguments as well as a functional 
renormalization of the vortex-vortex interaction energy to calculate the effective Gibbs free energy 
on the scale of the mean flux line distance. We discuss the resulting low field phase diagram and 
make quantitative predictions for pure BiSGCO (Bi2Sr2CaGu2 0s)- In the case with impurities, the 
Gibbs free energy is calculated on the basis of scaling arguments, allowing for a semi- quantitative 
discussion of the low-field, low-temperature phase diagram in the presence of impurities. 

PACS. 74.6O.Ec Mixed state, critical fields, and surface sheath - 74.6O.Ge Flux pinning, flux 
creep, and flux-line lattice dynamics - 74-72. Hs Bi-based cuprates 



1 introduction i/c2j respectively, because of the softening of cge- This ap- 

pUes in particular to high-Tc materials with their elevated 



I 

O . Conventional type-II superconductors show in addition to transition temperatures and their pronounced layer struc- 

^ the flux repulsing Meissner state a second superconduct- ^^^^es ||. At present, it is not clear whether the transition 

> ' ing (Abrikosov) phase in which the magnetic induction ^ ^^^^ normal phase at high field happens m these mate- 

■r^ B enters the material in the form of quantized flux lines "^Is via one or two transitions. However, melting of the 
?S (FLs) which form a triangular lattice. Each FL carries lattice has clearly been observed experimentally 0. 

i a unit flux quantum = hc/2e. The Abrikosov lattice 
■ - - ' is characterized by a non-zero shear modulus cee, which At low fields a first order melting transition to a liquid 

vanishes at the upper and lower critical field, H^^ and phase and a change in the critical behavior of B has been 

Fci, where continuous transitions to the normal and the predicted some time ago by Nelson||. Quantitatively the 

Meissner state, respectively, occur. In his mean-field solu- influence of thermal fluctuations is described hy a thermal 



tion Abrikosov treats FLs as stiff rods. Close to the lower length scale Lt = <PII{1&tt'^T) k. 2cmK/r |§. Lt has a 

critical held H^^ , their interaction becomes exponentially simple physical meaning: an isolated flux line segment of 

weak and hence the FL density ^ qj^^ vanishes as length Lt shows a thermal mean square displacement of 

I ln/i|-2 where h = {H - H,,)/H,, denotes the reduced order of the London penetration length A. Besides a 

field strength [Q| shift of , large scale thermal fluctuations lead close to 

Thermal fluctuations roughen the FLs resulting in a ^o an entropic repulsion ~ (AVLt^)' between FLs 

possible melting of the Abrikosov lattice close to i/^ and which dominates over the bare interaction for small h and 

hence i? ~ /i Here and below, all FL interactions are 

^ Email address: av@thp.uiii-koeln.de measured in units of £0 — (^o/47rA)^ = Lj-T j}?. 
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More recently, Blatter and Geshkenbein g found that 
in anisotropic or layered superconductors short scale fluc- 
tuations give rise also to an attractive van der Waals (vdW) 
interaction 0. For FLs separated by a distance R the 
strength of this interaction is of the order —X^/{dLTR'^) 
for A < i? < d/e and of the order -\*^ /{sLtR^) for 
d/e < R < X/e. = m/M ^ 1 denotes the anisotropy of 
the material with m and M the effective masses parallel 
and perpendicular to the ab plane, and d the interlayer 
spacing. A and A/e are then the screening lengths parallel 
and perpendicular to the layers, respectively. 

The competition among the bare, the entropic and the 
vdW interactions leads to an interesting phase diagram at 
low B values. In particular. Blatter and Geshkenbein 
find at low T a first order transition between the Meissner 
and the Abrikosov phase. 

So far fluctuation effects have been discussed for a 
clean superconductor. It is well known, however, that in 
type-II superconductors FLs must be pinned in order to 
prevent dissipation from their motion under the influence 
of an external current. Therefore, besides the thermal fluc- 
tuations one has to take into account the effect of disorder. 
Randomly distributed pinning centers lead indeed to a de- 
struction of the Abrikosov lattice but as has been re- 
cently shown, for not too strong disorder FLs may form a 
(Bragg-) glass phase which is characterized by quasi long- 
range order of the FL lattice |9|. For low B this phase 
undergoes a melting transition to a pinned liquid state. 
Inside this phase, disorder induced effects are expected to 
dominate over those of thermal fluctuations for sufficiently 
low T. The influence of the disorder fluctuation induced 
forces between the FLs and the consequences for the low 
B phase diagram - which deviates substantially from that 
found in Ref. ^ for pure systems - have been considered 
in @. 

In this paper, we address several issues. We start with a 
review of the interactions between FLs and the derivation 
of the properties of an isolated FL in section ^. Correlation 
functions, both in the thermal and in the disordered case, 
are calculated. These are needed for the calculation of the 
van der Waals interaction in section ^. There, we first de- 
velop an intuitive picture to derive the vdW energy driven 
by thermal fiuctuations and by disorder. In the presence 
of strong disorder, detailed level statistics of the random 
impurity distribution are necessary. In the second part of 
this section, a thorough derivation of the van der Waals 
interaction by a systematic expansion of the free energy 
of two fluctuating FLs is presented. 

The consequences for the low field phase diagram of 
anisotropic high- Tc superconductors are considered in sec- 
tion ^. In the thermal functional renormalization 
group is used to calculate the effective Gibbs free energy 
on the scale of the mean distance £ between FLs. In par- 
ticular, the bare interaction between flux lines, which is 
given by the superposition of the bare repulsion and the 
vdW attraction, is renormalized by integrating out ther- 
mal fluctuations on scales between A and £. The results 
of this calculation are compared to an expression for the 
Gibbs free energy based on scaling arguments. In this lat- 



ter approach, the contribution from the vdW interaction 
can be estimated only up a to numerical factor determin- 
ing its amplitude; this factor will be quantified using the 
results from the renormalization procedure. In the case 
with disorder, we give an expression for the Gibbs free en- 
ergy that is based on scaling arguments, and which allows 
for a semi-quantitative discussion for the low-field, low- 
temperature phase diagram in the presence of impurities. 



2 Single chain properties 

In this section we discuss properties of an isolated flux 
line (FL). Starting with a vortex lattice picture with the 
general form of the interaction between the FLs, we con- 
sider the limit of a very dilute lattice where a single chain 
approximation is good enough. A single FL in the dilute 
limit is well described by a dispersive stiffness constant 
ei{kz), the origin and different limiting properties of which 
are discussed. Deriving the Hamiltonian for a FL in an 
isotropic material, we generalize it for a FL in the pres- 
ence of impurities with its full dispersive stiffness constant. 
Finally, we discuss the properties of several correlation 
functions in the presence of impurities. These correlation 
functions will be used in the derivation of the van der 
Waals interaction in section 0. 



2.1 Interactions and elastic constants 

In reality the flux lines are not straight; they are distorted 
from their equilibrium positions. The displacement with 
respect to the equilibrium position = (a.i,z), where i 
denotes a site on the planar lattice, is represented by a two 
dimensional vector Ui{z) In the low field limit, where 
the London theory is applicable, the interaction energy 
between N vortices is in general given by 




where = + Ui{z) is the position vector of the dis- 
torted vortex segment, (a/?) S {x,y,z) indicate differ- 
ent components. In (|^) the terms with i = j have been 
included. These terms correspond to the interaction be- 
tween the segments of the same vortex line and there- 
fore contribute to the self energy of the line. In isotropic 
materials the interaction is exponentially damped and is 
of the form e~'^^^/r with a short distance cutoff at the 
coherence length r — ^. The high temperature oxide su- 
perconductors are essentially layered with the conducting 
CuO planes parallel to the ab plane. A large anisotropy 
between the c axis and the ab plane due to this layered 
structure leads to many unusual properties compared to 
the isotropic materials. 

The tensorial London potential between the vortices 
is conventionally expressed through its Fourier transform 
defined as Vif{r) = AnX'^ J d^k/{2TT)^ e'^'^V^fik) In 
our case, however, it turns out to be more convenient to 
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express it in terms of its partial Fourier transform with 
respect to the z-direction only, as defined by 



T/int 



(R, 



1 



47rA2 



dz e 



(2) 



where r = (R, z). This choice has the advantage of an ex- 
plicit dependence on the distance R between FLs. Choos- 
ing R = (i?, 0) and i? ^ A, the only contributions to the 
tensorial London potential that do not vanish by symme- 
try are 



1 



27rA2 



Xo(yTTA2fcfi?/A) (3) 



and 



int 
VV 

1 £ 



(4) 



'2TrRX^l + X^k 

Kn{x) is the modified Bessel function of the second kind 
of v}^ order. In the extreme anisotropy limit e — > 0, this 
latter equation simplifies to 



1 



1 



(5) 



In (|^) to (^), a short distance cutoff at the coherence 
length ^ is implied. In (^) and (||) subdoniinant terms 
of order e~^^-^ have been neglected. 

Different elastic moduli of the vortex lattice can be 
obtained from by expanding in small displacement u. 
For very small field H < H^^ , the shear and compres- 
sion moduli cgg and cii(fc) decay exponentially with the 
lattice spacing £ and can be neglected in a very dilute 
limit. The tilt modulus C44(fc) is only weakly dispersive in 
isotropic materials ||l^. In anisotropic materials, though, 
it has a much stronger dispersivity [p"3[|^ which results, in 
the single FL limit, in a strongly dispersive elastic stiffness 
constant £i{kz) — l?'c%^{kz) consisting of two parts. 



■In 



1 + X^kt 



■Hl + klX^), (6) 



2A2A:; 



where Kc = i^/£ with the Ginzburg-Landau parameter 
K = A/^. c%j^{kz) is the part of the tilt modulus that orig- 
inates from the self-energy of the vortices. The first part 
in (|6[) represents the Josephson coupling between the lay- 
ers and is absent in the limit e ^ 0. The second part 
describes the contribution from the electromagnetic cou- 
pling, which remains finite in this limit. In the extreme 
short wavelength regime Xkz > 1/e the contribution due 
to the Josephson coupling dominates over the electromag- 
netic part and the line stiffness is given by the first term in 
(^. In the intermediate regime 1/X < k^ < 1/eA, the elec- 
tromagnetic interaction dominates, and the elastic con- 
stant is given by the second term in (^, hence the stiffness 
is highly dispersive in this regime. In the long wavelength 
limit, the line stiffness reaches a constant. 



£i{kz 



' In ni 



for Xkz < 1. 



(7) 



In order to obtain an estimate for the contributions 
from the Josephson and the electromagnetic coupling to 
the stiffness constant at small length scales we use kz = 
Tr/d, where the layer spacing d is the lowest relevant length 
scale. Using the parameters A = 2000 A, ^ = 20A, d = 
15 A and e = 1/300 |], suitable for BiSCCO (BiaSra- 
CaCuaOs), we find that the contributions from the elec- 
tromagnetic and the Josephson couplings are about equal 
in magnitude at this wavelength, while for lower kz, the 
electromagnetic part dominates. Hence, for BiSCCO the 
Josephson coupling may be entirely neglected. For mate- 
rials with a larger value of e, on the other hand, kz = 7r/d 
belongs to the extreme short wavelength limit [Xkz > 
1/e) where the weakly dispersive Josephson coupling dom- 
inates. 



2.2 Hamiltonian for a single Flux line 

The line stiffness £i{kz) obtained above can be used to de- 
scribe FLs as elastic strings with a strongly dispersive elas- 
tic constant. In practice it is difficult to predict the prop- 
erties of the FL at all wavelengths. Often in the following 
we therefore aim at large or small wavelength features de- 
pending on our interest in distortions on scales larger or 
smaller compared to A. In the extreme long wavelength 
case where FLs are elastic strings with a dispersionless 
stiffness constant, the underlying lattice structure can be 
ignored. In the opposite limit the layered structure and 
the anisotropy become important Depending on the 
ratio Tc = 2Q/d'^, where — £(, is the coherence length 
in the direction parallel to the c axis with e denoting the 
anisotropy parameter, one can distinguish two cases, the 
continuous anisotropic limit and the extreme decoupled 
discrete limit. For Tc ^ 1 the coherence length is large 
enough to neglect the discrete layered structure. In that 
case the continuous anisotropic Ginzburg-Landau descrip- 
tion is applicable. Tc ^ 1 is the decoupled limit where the 
Lawrence-Doniach model is more appropriate. In the 
latter case, instead of the continuous description, the vor- 
tex is often viewed as a stack of pancake vortices residing 
in layers but interconnected by Josephson strings between 
two successive CuO planes, as sketched in Fig. ^ The ef- 
fect of the thermal energy or of random impurities is to 
displace the pancakes from their aligned position. Due to 
the dispersion these local displacements and their corre- 
lations will be different from the long wavelength distor- 
tions. It is however important to point out that the contin- 
uous anisotropic Ginzburg-Landau theory often provides 
a good description of the layered structure jl5| unless spe- 
cific thermodynamic properties, where the 2d structure is 
more important than the 3d bulk material, are investi- 
gated. 

Taking into account the dispersivity of the stiffness 
constant of the FL, the Hamiltonian for a single FL in the 
presence of random impurities can be defined as 



Hi{n}) = 



£i{k) 



k 



PukU-k+ dzepin(u{z),z).{8) 
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Fig. 1. A single vortex line, carrying the flux $0, in a strongly 
layered superconductor. The flux line consists of an array of 
pancake vortices threading the individual superconducting lay- 
ers which are separated by a distance d; the pancakes are in- 
terconnected by Josephson strings. 



Here only the bending energy part has been expressed 
in momentum space with the Fourier transform of u(z) 
defined as 

VI. 



Ufc = 



dz 



^u(z). 



(9) 



Here and below, k denotes the fcz-component of the full 
wave vector k. 

The random pinning potential is assumed to be Gaus- 
sian distributed with epin(u) = (a nonzero mean value 
Epin simply shifts the Hamiltonian by a constant and is 
hence inconsequential) and 



epi„(u, z)£pin(O,0) 



rp3 



(10) 



with k{x) = 1 for x ^ 1 and k{x) ~ (l/a;^)lna; for x ^ 
1, respectively Q|. The overbar denotes the average over 
the disorder distribution. The displacement caused by the 
quenched disorder is strongly suppressed due to thermal 
smoothening for T > Tdis. The characteristic energy Tdis 
can be related to the local shift in the critical temperature 
Tc which is induced by an oxygen impurity density Ui by 



7^3 



0.4 



Tr. drii 



1 



1 



(11) 



where t = T/Tc This result holds in the case of non- 
optimal doping where dTc/ dui ^ 0; if this derivative van- 
ishes, Tdis can be similarly related to the second derivative 
d^Tjdn-i p. 

Statistical mechanics of a single FL in a random poten- 
tial shows that the disorder is always relevant for d < 2. 
This is reflected through an enhanced averaged FL wan- 
dering on large scales L > Ldis 



([u(i)-u(0)]2) ^ A(L/LdiO^ 



(12) 



with C>l/2for(i<5 [Q. Here and in the following 
the angular bracket indicates thermal average. Several ex- 
act treatments at d = 2 lead to C, — 2/3. Disorder is 
marginally relevant at d = 3. Numerical simulations and 
recent analytical calculations show C = 5/8 ||l^, 18|. In the 
presence of disorder vortex segments of length Ldis are 
independently pinned. At T = 0, this can be seen by min- 
imizing the energy consisting of elastic and pinning con- 
tributions. Here the pinning energy is estimated within 
perturbation theory. This provides 



^dis — £oA K 



2^1/C-2 



/Tdis. 



(13) 



For T < Tdis, the thermal fluctuations can be ignored and 
the above estimate of Tdis is still valid. For T > Tdis, on 
the other hand. 



LdisiT) 



T 



■ exp 



(T/Tdis)' 



(14) 



shows an exponential growth with the temperature as a 
consequence of the marginal relevance of the disorder at 
3 space dimension p9[ |. In the pure thermal case ^ = 1/2, 
and Tdis is replaced by [|j 



(15) 



2.3 Correlation functions 



With the ideas sketched above we now study several cor- 
relation functions in the presence of quenched impurities. 
Whereas in a pure system (ufc) — and hence (ufcU_fc) is 
the appropriate correlation function, in a random system 
one has to distinguish between the connected correlation 
function 

CT(fc) 



and the disconnected correlation function 



Cdis(fc) = (Ufc)(u_fe). 



(16) 



(17) 



We first study Cdis(fc) and near the end of this section we 
show that the disorder does not affect Cxik) obtained in 
the pure case. 



2.3.1 Short wavelength limit of Cdis(fc) 

We first consider the short wavelength behavior of Cdis(fc). 
In this limit the problem essentially simplifies to obtaining 
the typical displacement of a single pancake. A perturba- 
tive technique and an Imry-Ma type argument are applied 
to estimate the single pancake di spla cement. 

As mentioned above in section 2.1 , for highly anisotropic 
superconductors with very small values of e ^ 1 (such as 
BiSCCO), the highest relevant wave vector Tr/d is smaller 
than 1/Xe, with the consequence that in the short wave- 
length limit Afc ^ 1 it is sufficient to consider the disper- 
sive elastic constant e;(fc) = (so/A^fc^) In(Afc). This term 
represents the contribution from the electromagnetic cou- 
pling in the full expression m) for the line stiffness. A 



Andreas Volmer et al.: Van der Waals interaction between fiux lines in High- Tc Superconductors 



5 



closer look at the elastic energy term in (||) with the dis- 
persive elastic constant reveals that, apart from a logarith- 
mic factor, the pancakes on different layers are essentially 
uncoupled. Thus rewriting the correlation function in real 
space 

(Ufe)(u_fc) = ^1^^ d^'e-''=(^-^')(u(z))(u(z')), (18) 

we utilize the above assumption of independent pancake 
displacement, i.e. 



{u{z)){u{z'))=uUiz-z'), 



(19) 



where 



/(^) 



1 for \z\ < d/2 

for d/2 < \z\ < A. 



This leads to 



Cdis(fc) ~ (^2 ) d for K Afc < nX/d. 



(20) 



The mean square displacement of a pancake in the pres- 
ence of quenched impurities can be obtained from the 
Hamiltonian (^ which is, in the short wavelength limit, 
a superposition of the following two parts: (i) a parabolic 
elastic potential, Ec\{u) « {v?d/2) A:2e/(A:)| j,^^/^, which in 
the case considered here (where the electromagnetic cou- 
pling dominates) reads 



Eo\{u) = i^ln(7rA/d)|j, 



(21) 



and (ii) a contribution from the disorder, averaged over 
the pancake height rf. 



£pin(u) = dz epin(r -I- u(z), z) 
Jo 



(22) 



where r is the equilibrium position of the pancake. 

It is useful to re-express this Hamiltonian in dimen- 
sionless variables t = z/d and u = u/^. In terms of these 
dimensionless variables. 



Ti Sod 



2 ln(7rA/d) X 



(23) 



^ + A"\n\/d) j\t epi„(u,t)| , 



where epin(u,f) is the pinning potential scaled to have a 
variance equal to unity, and we have defined 



A{x) = Z\o 



ln^(l-|-a;2) 



^0 



V eoA 



(24) 



From ( p3D one concludes that in general, at finite temper- 
ature and finite disorder strength, 



-'pv) ~ 



A{n\ld), 



Td 



eeiin/d) 



(25) 



where and ^ is a function of dimensionless variables for 
the disorder and the temperature. 

In the absence of disorder {Aq = 0) fluctuations are 
driven by thermal energy and we have 



("pv)th 

e 



2T 



Sod ln(7rA/(i) 



(26) 



At r = and finite disorder (T^is > 0), the FL fluctuates 
to take advantage of the impurities, and the ground state 
displacement in this case is 

ulje - F{A{^\ld)) EE <Z>(Z\(^A/d), 0). (27) 

For very low temperatures T <C Tdis, we can approximate 
the mean square displacement as 



*pv/ 



u2 =e2F(Z\(7rA/d)). 



(28) 



In order to determine the explicit functional form of F{A), 
we take a recourse to a perturbative approach valid for 
small displacements w <C ^. In order to go beyond the 
perturbative regime, an Imry-Ma type scaling argument 
and the result from numerical simulations are discussed. 
In general we expect a power law form of 



F{A) - zi" 



(29) 



and we attempt to get the value of rj. 

In the perturbative scheme we proceed with the re- 
duced Hamiltonian (^3|) for a single pancake vortex. Ex- 
pansion of the pinning energy in small displacement u 



epin (U, t) = Spin (0, t) + U • VSpin (u, i) | 



(30) 



and minimization of the total energy with respect to the 
displacement leads to the force equation 



A^'\TT\/d) / dtV£pin(u,i)| 



(31) 



Averaging over the dimensionless distribution £pin(u, i) 
directly leads to the disorder induced short length scale 
displacement 

^,^eA{T:\/d). (32) 

Thus, within the perturbative approach -q — 1. 

Due to the narrow range of applicability of the pertur- 
bation technique, a more general nonperturbative scheme 
is needed for A{T:X/d) 3> 1. This can be done by an Imry- 
Ma or variational approach ||2^ , which briefly proceeds as 
follows. Consider the dimensionless part (in curly brack- 
ets) of the single pancake Hamiltonian (p3|). The displace- 
ment of a single pancake by an amount u is associated 
with the elastic energy cost u2/2. This energy cost has to 
be compared with the possible gain in pinning energy cor- 
responding to the adjustment of the pancake position to 
a minimum of the random potential within the area ttv? . 
This energy is obtained from the second term in (^) and 
is of order A^/'^{'k\/ d), where we have used that £(u, t) 
has been scaled to unity, i.e. 



(itepi„(u,t) / rft'epin(u',i') = 1. 



(33) 
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Balancing the two energies we find 



2.3.3 Thermal fluctuations Ct 



(34) 



corresponding to the exponent rj = 1/2. 

In the derivation above we have neglected the fact that 
the gain in pinning energy is not independent of u. With 
e(u,f) obeying a Gaussian distribution, it is in fact pro- 
portional to ln^^^(u^/2y^) Taking account of this 
dependency leads to a confluent logarithmic term of the 
form F{A) - zi^/^ ln"^/^(Z\). Nevertheless, we will only 
use the leading scaling behavior of F{A) in the rest of the 
paper. 

A numerical estimate for BiSCCO with the parame- 
ters as used before and Tdis = 45 K yields A{-K\/d) ~ 300, 
hence the Imry-Ma result should apply, predicting a typ- 
ical displacement Upv/C^ of order 20. 

Finally, we mention the result from a numerical study. 
We have simulated the energy on a 2-d lattice by super- 
posing the elastic energy quadratic in the position r on the 
lattice and a random, Gaussian distributed energy V{r) 
pO| . A log- log plot of the mean square position of the min- 
imum energy site, obtained after a large number of sam- 
ple averaging, versus the width of the distribution, shows 
roughly a straight line with a slope that corresponds to an 
exponent 77 « 4/9. While this may not reflect the correct 
asymptotic scaling behavior, it gives a good representation 
of the scaling in the regime of physically accessible values 
for Z\, and is consistent with the confluent logarithmic 
term mentioned above. 

In conclusion, we find in the short wavelength and low 
temperature limit 

Cdis(fc) = w^^d for l<Afc<7rA/d, (35) 

where ^ w C^F(zi(7rA/d)) and F{A) ~ Z\'' with 77 = 1 
for Z\ < 1 and r/ w 1/2 for Z\ > 1. 



2.3.2 Long wavelength limit of Cdis 

In the long wavelength limit Afc <C 1, we are aided by the 
known scaling form of the mean square displacement ( p^ 
and the value of Q from numerical simulations at d = 3. We 
proceed with an Imry-Ma argument as described above, 
with the replacement of irX/d by Afc. For Afc <C 1 we have 



Zi(Afc < 1) ~ Ao/{\kf 



Therefore using 



l/L 



dk 

2^ 



Cdis(fc) -e' 



l/L 



dk 

27rfc 



(36) 



F{A{\k)) (37) 



together with (|T|), we find F{A{>Jt)) - [A{\k)Y^/^. This 
leads to Ldis in agreement with (n3[) . The above treatment 
is valid for very low temperature T <C Tdis- At higher tem- 
peratures the crossover length scale grows exponentially as 
given by Eq. (|lj) . 



The other quantity of interest is the disorder averaged 
second cumulant of thermal fluctuation Cxik) as defined 
in ([l7|). In the pure system the single vortex fluctuation 
amphtudc is (ufcU_fc) — T/ei{k)k'^. In the presence of dis- 
order, Cxik) remains unaltered, i.e. 



Crik) = 



ei{k)k^ 



(38) 



This can be readily seen from the discretised single vortex 
Hamiltonian in a random potential by adding a source 
term XiUi. Derivatives of the partition function with 
respect to Ai lead to various correlations. The above result 
follows by rewriting the Hamiltonian with a new variable 
Uk — TXk/eik'^ in terms of which the random potential is 
assumed to have the same distribution 



3 Van der Waals attraction 

In the first part of this section we discuss the origin of 
the van der Waals (vdW) attraction considering a sim- 
ple physical picture of distortion of pancakes and sub- 
sequent formation of pancake-anti pancake dipole. This 
simple approach correctly accounts for the vdW attrac- 
tion in the pure system. Using the same picture, we also 
discuss the vdW attraction in the presence of impurities in 
the medium - a problem more complicated than the pure 
case due to the averaging over the impurity distribution 
which has to account for the pres ence of metastable states. 
In the second part (section 3.2), we derive the vdW at- 
traction using the statistical mechanics of interacting flux 
lines. The power law behavior of the vdW interaction as 
a function of the FL distance R differs for the Josephson 
coupled and uncoupled cases due to the very nature of the 
interaction. 

For a stack of only electromagnetically coupled pan- 
cake vortices the interaction between two pancake vor- 
tices in the same layer separated by a distance R > ^ 
is given by V{R) = 2c?eo ln(i?/^). Pancakes on different 
layers attract each other with a strength reduced by a 
factor d/X. Josephson coupling between the layers leads 
to an additional contribution to the free energy that de- 
pends on the phase difference of the superconducting order 
parameter on two successive layers. A comparison of the 
Josephson and the electromagnetic interaction energy 
shows that for separations R < d/e the latter is much 
stronger, whereas the Josephson interaction dominates in 
the regime R > d/e. Therefore the interaction between 
two Josephson coupled pancakes in the same layer is given 

by 



V{R) = 



2deo In — , 



sR d \ 
T ~ leRj 



(39) 



2deo ( ^ - T-^ \ ,d/e<R< X/e. 
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For two straight stacks of pancakes placed at a distance 
i? ^ A, the mutual attraction from different layers and re- 
pulsion from the same layer compensate each other leav- 
ing a residual interaction ^ exp(— _R/A) as it is in the 
Abrikosov lattice. 



There, however, we will show that a different argument, 
using level statistics, leads to an expression for the vdW 
energy that is basically identical to (^3|). Hence, this lat- 
ter expression fully applies both in the thermal and in the 
disorder dominated regime, justifying the thermal anrf dis- 
order average over which has been anticipated here. 



3.1 An intuitive picture 

A FL with a single pancake distorted due to thermal en- 
ergy or disorder can be visualized as a straight FL with 
a dipole constructed out of the displaced pancake and an 
anti-pancake at the vacant position left behind by the dis- 
placed pancake. The dipole thus formed induces another 
dipole on the neighboring vortex by virtue of the inter- 
action between the pancakes of different vortices and a 
finite restoring force of the same stack. These fluctuation 
induced dipoles interacting with the long range interac- 
tion originating from ( |39| ) lead to a long range interac- 
tion between vortices very similar to the quantum fluc- 
tuation induced attraction in atomic system. In analogy 
with the latter system these fluctuation induced attraction 
are called van der Waals interaction. 

If two pancakes in the same layer of stacks 1 and 2 
are displaced by an amount ui and 1*2, respectively, the 
dipole-dipole interaction energy between them is given by 

f/i2 = ^[ui-U2-2(ui-n)(u2-n)]. (40) 

Here n denotes the unit vector along the line connecting 
the two vortices. The force exerted on the pancake of stack 
2 due to a dipole on stack 1 (in the same layer) is 



2£orf 
i?2 



Ui — 2(ui • n)n]. 



(41) 



The displaced pancake on stack 2 on the other hand ex- 
periences an elastic force due to all the other pancakes in 
different layers. The restoring force £22 = — Vu2£'oi(u2), 
where £'ci(u2) is the parabolic elastic potential ( |2l] ) act- 
ing on the pancake vortex, is 



-EoTiT In {nX/d) U2. 



(42) 



In the absence of impurities, the position of pancake 2 is 
(in equilibrium) exactly determined by the force balance 
fi2 + f22 = 0, which leads to U2 = -(2AV-R^ ln(7rA/d))_x 
[ui — 2(ui • n)n]. Substituting this result back into ( pO| ) 
gives the van der Waals interaction 



4£o 



A^ 



ln(7rA/d) i?4 



(43) 



Below, we will first evaluate this expression in the case of 
thermal fiuctuations only, and then reconsider the force 
balance argument in the presence of impurities. For weak 
disorder, we will find that this argument is still valid be- 
cause the random forces are much weaker than the restor- 
ing force while it fails in the strong disorder regime. 



3.1.1 Pure thermal case 

For the thermal case in pure systems, we find the vdW 
attraction by plugging the result ( |2^ ) for thermally driven 
fiuctuations into Eq. (E^), which leads to 



4eo 



T (\ 



vdw 



In (TrA/d) <i£o \R 



(44) 



This, apart from a numerical factor, is in agreement with 
the vdW attraction for pure superconductors derived in 

i- 



3.1.2 Disorder case at T = 

In the case of disorder-induced fluctuations, a more sub- 
tle treatment is necessary. First, let us distinguish be- 
tween the weak and the strong disorder limit, defined by 
A = AinX/d) ^ 1 and Z\ ^ 1, respectively. For conve- 
nience, we define the total pinning energy £^^l^{\ii) acting 
on pancake vortex i as in (p^). 

In the weak disorder regime {A <ti 1), the pancake 
position is determined by the force equation (|3l|), lead- 
ing to ensemble fiuctuations u^^ ~ ^'^A (^^ which are 
hence much smaller than the disorder correlation length 
^. Furthermore, the restoring forces exerted by the elas- 
tic coupling to the other pancakes in the same vortex (^) 
are much stronger than the random forces Vepjj^(ui). This 
can be seen by estimating the ratio between the second 
derivative of the random energy and the elastic constant 
£odln(7rA/(i)/A^ in a way analogous to (BSh, resulting in 



A2 



dso ln(7rA/(i) 



pm 



z\ < 1. 



(45) 



Hence, the force balance argument leading to the estimate 
( p3| ) fully applies in this regime. The van der Waals energy 

is obtained by substituting the perturbative result ~ 
^^Z\ for the correlation function into (^), yielding 



{R) 



-So 



k2 IninX/d) 



for Z\ < 1. (46) 



In the opposite limit of strong disorder {A ^ 1), the 
fiuctuations are of order ~ ^^^1/2 (jgj.ivg(j from the 
Imry-Ma argument (p^). In this limit, a force equation of 
the type used above is not sufficient to determine the re- 
sponse of the pancake position to dipole-dipole forces due 
to the presence of a large number of metastable states. 
Instead, we have to consider the global energy minimum. 
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In particular, in the absence of the dipole-dipole interac- 
tion Ui2 ( ^0|) the equilibrium position of each pancake is 
determined by the position of the global minimum of the 
superposition 

E,i{u,)+e^;}^{u,), (47) 

with Eei{ui) as defined in (|2l]). This minimum will be lo- 
cated at a typical distance of order ^A^^^ from the vortex 
center at u.^ = 0. Adding now the dipole-dipole energy 
Ui2 to (^), we look for the positional change 6u2 of pan- 
cake 2 induced by the interaction with pancake 1, keeping 
Ui fixed. The change 5u2 will lead to a modified dipole- 
dipole interaction U12 U12 +5Ui2, from which the vdW 
interaction can be derived as V^^'^ — SUi2/d, where the 
overline again denotes averaging over the disorder. 

For the single pancake in stack 2 we are focusing on, 
two different scenarios can emerge. First, for large U12, the 
original position U2 may become metastable if the energy 
landscape is tilted strongly enough (note that U12 is linear 
in U2) in order that another, formerly local, minimum be- 
comes the new global minimum; this scenario is illustrated 
in Fig. ^. If, on the other hand, U12 is so weak that the 
global minimum remains at the same pinning center, one 
can again apply perturbation theory to find an estimate 
for 5u2. 

First, let us have a closer look at the weak U12 case. 
Being in the strong disorder limit (A ^ 1), the esti- 
mate (^) shows that the restoring forces originating from 
the pinning center are now much larger than the elas- 
tic forces. Hence, an additional force f applied to pan- 
cake 2 will lead to a response U2 ^ U2 -I- Su2 with w 
—dadi3e'p^^{u2) ■ Su2,p , neglecting the elastic forces. Substi- 
tuting the dipole-dipole force (|4^) for f , we find a typical 
displacement Su2 (A/i?)2[l/Z\i/2 ln(7rA/d)][^]i/2 which 
is much smaller than ^ for i? > A, and hence a typical en- 
ergy change 6U^2^^ of order 

Sod — ,A^ 



i?2 



Ui ■ 6U2 



-Sod- 



1/2 _ 

?• (48) 



i?4 ln(7rA/d) 

We have neglected orientational dependences here, focus- 
ing on amplitudes. Notice that the factor A^^^^ compen- 
sates for the disorder scaling of the fiuctuation square 

« ^2z\^/2. Hence, if there were no pancakes in the 
stack whose optimal position changes under the influence 
of U12 (below we will find that this assumption is actually 
wrong), this would lead to a weak, disorder independent 
van der Waals interaction. 

The second possible scenario emerges when U12 is strong 
enough to induce a change of the order of the deepest en- 
ergy minima for the pancake in stack 2. In this case, it 
will induce a shift Su2 of order i^zi^/'* ^ ^, since the new 
minimmn may be located anywhere in a distance {u^Y^^ 
from the center. Hence, 



(49) 



where U12 is taken to be a typical dipole-dipole energy at 
the distance R, 

2eod 




Fig. 2. In the upper half, a one dimensional cut through a 
typical energy landscape, as given by (^), is shown; the bare 
elastic energy is shown as a dotted line. The two deepest min- 
ima happen to be located on opposite directions from the ori- 
gin, in a typical distance of order ^^'^ — C'^^^''- In the lower 
half, a linear energy contribution Ui2{u2) is added, sketched as 
a dashed-dotted line. Note that the center of the new parabola 
Eei{u2) + Ui2{u2) has been shifted to the left hand site. The 
addition of U\2 changes the order of the minima, hence this sce- 
nario corresponds to the 'strong [/12' case; the deepest minima 
are marked by small circles. 



So far, we have focused on a single pair of pancake 
vortices in the same layer. Now, a vortex line consists 
of a whole stack of pancake vortices. In a certain frac- 
tion of these, the strong U12 scenario will apply, where 
the dipole-dipole interaction induces a positional change 
5u2 ^ ^; the rest will stay within its old minimal energy 
position, giving only a marginal contribution to the vdW 
interaction. By estimating this fraction r^/\ via level statis- 
tics of the Gaussian random distribution, we can derive an 
estimate for the strength of the vdW energy as 



* vdw ^ 



\ (vA5Uf^°''^ -f (1 - ry4)5[/-<=-k^ 



(51) 



c/i2 = ;7i2(|u2|«e^i/4 



i?2 



(50) 



Below, we will find that the second term in (51), emerging 
from those pancakes with 5u2 <C ^, can be neglected. 

We will now calculate the fraction rj^ of 'active' pan- 
cakes in the stack. First, consider once again the situation 
in the absence of dipole-dipole forces, where the two pan- 
cakes take their optimal positions Ui and U2. The dipole- 
dipole interaction will induce a force ~{2£od/ E?)h on 
pancake U2, with b = Ui — 2(n-Ui)n (note that = u^). 
Let us divide the area of size ^2/\i/2 which is accessible to 
U2 in two half circles, the dividing line being perpendic- 
ular to b. In half of the cases, the deepest minimum will 
be located {before the addition of the dipole-dipole force) 
in the half space where the energy landscape is shifted 
upwards, as shown in Fig. ^. These pancakes can proba- 
bly profit from U12 by taking a new optimal position in 
the other half space. With the disorder correlation length 
given by ^, the energy landscape that is accessible to the 
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pancake can be represented by 



(52) 



effectively independent disorder realizations. We can hence 
represent the energy (^) in each half circle by a set of N/2 
Gaussian random numbers ej^ (? = 1,2 for the two sets, 
m = 1, . . . , N/2), each measuring the minimum of the to- 
tal energy in a cell of size f^. 

In Appendix H, the probability density pi2(fc) for the 
difference Se = st/- — el-, between the smallest values of 



the two sets, respectively, is derived; cf. Eq. (103). Given 
this probability density, the fraction of 'active' pancakes 
can be estimated as 



Ul2 



d{Se) pi2{6e) 



(53) 



This integral gives the probability that the energy differ- 
ence between the deepest minima of the two sets is positive 
and smaller than the typical energy gain Ui2. 

For small enough values U12, corresponding to dis- 
tances R between the two vortices larger than some dis- 
tance Rc which we quantify below, we can approximate the 
probability density pi2((5£) in the integrand of Eq. ( |5^ ) by 
its value at = 0, 



Pl2{0) U12 



lni/2(iV/2) U^2 

e 



2tt 



4^, (54) 



where we have used (104), and e is defined as 



Now, plugging the results from Eqs. (54) and 
(|l|), we find 



(55) 
f) into 



Kdw(-f^ ^ ^c) — —£o 



4 Z\i/2 ini/2(^i/2) 

R 



27r k2 ln(7rA/d) 



for Z\ > 1. (56) 



It is easy to check that the second contribution in (^ is 
smaller than the first one by a factor of and can hence 
be neglected. Re-substituting a factor of ^^Z\^/^ by uf in 
(56|), we have finally arrived at an expression for the van 
der Waals energy which is identical to ( p3| ) up to a factor of 
0(1) and a subdominant logarithmic contribution in the 
disorder strength A. It is consistent with th e re sult from 
the systematic derivation obtained in section |3.2| , again up 
to a numerical factor and the ln(Z\^/^) term. 

Finally, we note that for distances R smaller than Rc, 
rjA —>■ ^, leading, via (|50|), to the saturation value 



So 



(57) 



At distances R < Rc, the vdW attraction is however dom- 
inated by the direct repulsion. The expressions (^6|) and 



(|j match at R 
given by 



Rc, where the crossover distance is 



R^ 



41ni/2 (/i/47r) 



(58) 



ln(7rA/d) 

Rc has a simple meaning: For R Rc, the typical dipole- 
dipole energy U12 is of the same order as the typical energy 
difference between the two deepest minima from the two 
sets. Since the prefactor on the right hand side of ( |5^ ) 
depends only logarithmically on A and X/d, Rc will be 
of order A, so that expression ( ^6|) correctly predicts the 
vdW interaction for all relevant values i? > A. 

With the parameters for the different length scales 
characteristic for BiSCCO as cited above, together with 
the disorder strength A = A{'K\/d) « 300, N = \fA is as 
low as « 20. With this value for A, we find Rc ~ cA with 
5 w 1.1. 



3.2 Statistical mechanical approach 

The vdW attraction can be derived using a more formal 
statistical mechanical approach where we consider the free 
energy of two FLs fluctuating due to thermal energy and 
impurities. In the absence of disorder we obtain the vdW 
attraction due to the temperature induced fluctuations of 
the lines. 

Our starting point is the free energy in (|l|) that con- 
tains a self energy part Tq and a mutual interaction part 
T\-[A consisting of terms with i 7^ j. To proceed further it 
is convenient to write the mutual interaction explicitly in 
terms of components: 



47r 



dk dz\dz2 e 



(59) 



V;'f(R-Kui(zi)-U2(z2),fc) 



t\a{zi)t2(3{z2) (R + Ul(2l) - U2(Z2), fc) 



where we have used ds^ = (t(2;), 1) dz with t(z) = 9zu(z). 
Here and below, each integral over fc implies a factor of 
l/27r. We split T\^t into a longitudinal part T\\ with the 
term V^z and a transverse part T with the term propor- 
tional to t\a{z\)t2(}^Z2) . The partition function which is a 
path integral corresponding to the weighted sum over all 
possible displacements is 

Z{R) = J I?[ui(z)] V[u2iz)] exp[-.F(i?)/T]. (60) 

In the presence of impurities, the disorder averaged free 
energy can be written as 



F{R) = -T\nZ{R) = -Tln(cxph.Fi„t(i?)/T])o, (61) 

where the average (...)o is taken with respect to the self 
energy part Tq. The effective interaction between the FLs 
follows immediately from the cumulant expansion 



LVc^{R) « (.Fint) - [{Tl,)^ - (•^int)^]/2T 



(62) 



Our aim at this point is to look for a long range vdW type 
attraction in T4ff. 
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3.2.1 First cumulant 

It is expected that the fluctuation induced vdW attraction 
should originate from the second cumulant. We first show 
that no such contribution arises from the first order term 
in ®. 

Starting with the longitudinal part 

(62 r 

JT,, = -0 / dkdzidz2 e'''^'-'-'^'^ X 
" in J 

yif(R + ui(zi) -U2(z2),fc) (63) 

a straightforward approach is to expand the interaction 
term V^f{R + Ui(zi) - U2(z2), k) in u. The 0(1) term 



^2 

in 



^ J dz2 (R,fc = 0) 



(64) 



with K,'"'(R, fc) as given by (H), leads to a short range 
potential decaying exponentially for R > X. The higher 
order terms lead to higher order derivatives in R, as can 
be seen by inspecting the O(m^) term 

g/d...,..2E^^Kf(R,^)e-^--)x 

a/3 

{(Ulaizi) - U2aiz2))iuif3{zi) - 1*2/3(2:2))), (65) 

and hence again are of short range in nature. 

The 0(1) term from a similar expansion of the expo- 
nential of the transverse part 



ik(zi —Z2) 



^ J dkdzidZ2 (tl„(zi)t2/3(22))'Ki^*(R,/c)e 

(66) 

requires averaging over the impurity distribution for the 
flux line displacements. We approximate the averaging as 
{tia{zi)t2i3{z2)) « (tiaizi)) {t2t}{z2)) ■ Such a decomposi- 
tion makes sense only when the flux lines are far apart so 
that the preferable impurity sites of line 1 do not affect 
the conflguration of the line 2. With this approximation 
it is now evident that the transverse part vanishes since 
{tia) = 0. The 0{u) term contribution vanishes for the 
same reason. O(u^) terms require evaluation of several 
correlation functions. As an example we pick up one term 
and show that it vanishes due to the symmetry of the inte- 
gral. In the same way, it can be shown that also the other 
terms vanish. Let us consider the term 



2- — / dk dzi dz2 

An 



||e;(R,^)e 



{tla{Zi)t2(){Z2)uix{zi)u2x{z2))- (67) 

As before we assume that the averaging over the impurity 
distribution for FL 1 and 2 is uncorrelated in nature. The 
above correlation function is therefore equivalent to 



{tla{zi)uix{zi)) {t2p{z2)u2x{Z2)) (68) 

= - j d{h}kik3 (C?-(fci,fe) + CS:(fci,A;2)) X 
(C^"(fc3, fc4) + 0^^(fc3, fc4))e^('=i+'=^)-ie^('=3+fc4)^2^ 



where we have used C^f (fci, k2) 



{uia{ki)){uif3{k2)) and 



m a 



C^'^{ki,k2) = {Uia{ki)uif3{k2)) - (Mla(fcl)) ("1/3(^2)) 

straightforward gen eralization of the corresponding quan- 
tities ( |l6| ) and ( [1^ ) defined in section |^. Due to trans- 

lational invariancc, both C^-^{ki, k2) and C^^ {ki,k2) are 
diagonal in their arguments, i.e., they include a factor of 
5{ki-\-k2). Using this property of the correlation functions, 
the integration over zi, Z2 and k in (p37|) leads to 



2^ 

An 



dkidk^ kik; 



0) X 



{C^%k^) + CS,l{k,)) (Cf (fc3) + (69) 

where we have switched back to the notation as in ( |T^ ) 
and (p7|). Now, both Cdis(fc) and Cxik) are symmetric 
in k, which leads to an integrand odd in ki and k^. The 
integration regime given by ki, = —n/d...n/d, the 
integral hence vanishes by symmetry. 



3.2.2 Second cumulant 

Next we go over to the second term in (^2[). The pure 
longitudinal part again leads to only a short range 
interaction which has the same effect as before of renor- 
malizing the short range potential between the flux lines. 
The cross term of the longitudinal and the transverse parts 
contains (^10^2,3) and therefore vanishes after disorder av- 
eraging if the flux lines are assumed to be far apart. We 
are left with the contribution from the pure transverse 
part [{J-'j^)o — {T±)q]/2T, which has to be analyzed in a 
more detail. Indeed, as we will show the contribution to 
the vdW attraction per unit length can be found as 



1 



vdw 



2TL 



(70) 



To lowest order in the displacements, J-± reads 

^^^'^Jdkk^ uio.{-k)u2p{k) Ki^*(R, k) , (71) 

with V^^*(R-, k) as deflned in (|) or (|) for e > or e ^ 0, 
respectively. 

For the evaluation of (^0|), it is now necessary to eval- 
uate the correlation function 



C{k) = {uia{-k)u2p{k)ui^{--k')u2s{k')) 

-{ui^{~k)u20{k)){ui^{-k')u2s{k')). (72) 

Using the above decomposition for disorder averaging, this 
can be rewritten as 



C{k) = C^^-k, -k')C^\k, k') + C^i-k, -k') X 



(73) 



The first term in ( 73 ) involving only Ct leads to the vdW 
attraction found in the pure system. This follows from the 
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fact discussed in section g that the disorder does not affect 
Ct- Substituting this term in (|7^), we find 



a/375 



(74) 



Using the fact that C^"' and V™^ are diagonal in com- 
ponents, and substituting V™^ from (||), we find in the 
extremely decoupled limit 



1 



1 



2T (47r)2 (27ri?2)2 

£0 T /A_ 
ln^(7rA/d) dso \R 



dk 



(1 + A2fc2)2 \k^ei{k) J 



(75) 



In the last step, we have approximated the integrand by 
its value at the upper integration limit ir/d where the con- 
tribution to the integral is the largest. 

In the continuous anisotropic case one has to con- 
sider the more general form of the potential in (Q) . Ki (x) 
has the asymptotic properties Ki{x ^ 1) ~ 1/x and 
Ki{x 3> 1) ~ e~^/y/x. Due to the exponential decay of 
Ki{x) for large arguments, the vdW attraction is expo- 
nentially suppressed for R > A/e. For smaller 7?, without 
much loss we may set the upper limit of the integral to 
1/eR and use Ki{x) ~ 1/x in the integrand. Again the 
contribution is maximum at the upper limit. Thus we find 
the thermal vdW attraction in the continuous anisotropic 
case 



T 



vdw 



ln'^(7rA/ei?) eAeo 



for - < i? < - . 

e £ 



For R < die, Ki[x) ~ 1/a; in the whole integration range, 
so that ( [Tq ) applies in this regime. 

As expected, the vdW attraction in a pure system van- 
ishes as T ^ 0. Expressions (|7^) and (|7^) for the ther- 
mally induced vdW attraction have been obtained previ- 
ously in Ref . |^ . 

Using the fact that the correlation functions are diag- 
onal in components, the most general form of the vdW 
attraction in the presence of the disorder and at T ^ is 



* vdw 



1 



Txjd 



4T (4^)2 



dk fc* 
+ 2CT(fc)Cdis(fc)) 



(77) 



In order to see purely the effect of the disorder, one has to 
study the the term with Cdis in this equation. In section 
(2.3), we have derived expressions for Cdis(fc) both in the 
limit AA: ^ 1 and Afc ^ 1; the exact behavior at AA; of 
order unity is however unclear. Nevertheless, provided the 
roughness exponent < 1, one finds that the dominant 
contribution to the fc-integral comes from large wave vec- 
tors k < n/d. Hence, the vdW attraction can be derived 
by substituting the expression for Cdis from (^ , and pro- 
ceeding as in the thermal case. It is found to have the same 



_R-dependence as in the thermal case both in the extreme 
decoupled and in the continuous case, respectively. We 
give the result in form of an ad-hoc interpolation formula 
which subsumes both limits R <C d/e and R ^ d/e, 



^^dw 



^0 



d 



d + eR 



In- 



ttA 



-1-2JJ 



eR, 

(78) 

with Aq and 77 as defined in ( p^ and (^9|), respectively. 
Like in the thermal case, ( [78| ) holds only for R < A/e; for 
larger R, the vdW attraction is exponentially suppressed. 
Note that in the decoupled limit e ^ 0, the result from 
the intuitive approach (^6|) is reproduced up to numerical 
factors. 

Finally, we note that higher order terms in the dis- 
placement u in the expansion of the interaction potential 
in (59|) introduce derivatives of V^^* by components of R 
The 0{u.^) term, for example, yields a finite con- 



tribution at T = and finite disorder, but vanishes ' 
or i?^'', respectively, and can hence be neglected. 



R- 



4 Phase Diagram 

In the low- field limit, where the magnetic induction B is 
of order a few Gauss, the attractive van der Waals interac- 
tion between vortex lines has consequences for the phase 
diagram of anisotropic superconductors. With the exter- 
nal magnetic field H fixed, one has to analyze the Gibbs 
free energy density 



G(f;i/,r,Tdis) =F(£;T,Tdis 



eo In K - 



(79) 



which has to be minimized with respect to the mean FL 
distance i. h ^ {H - H^J/H°^ is the reduced magnetic 
field, H^_^ = 47reolnK/«?o being the unrenormalized lower 
critical field, and F{i; T, Tdis) is the free energy density. £ 
is related to t he m agnetic flux via B — 'Po/£'^. 



In section 4.1 
thermal case T 



> 



we will restrict ourselves to the pure 
with vanishing disorder (Tdis — 0). 
We will give a short overview over approaches applied 
previously to the determination of the effective Gibbs free 
energy den sity, a nd review a scaling approach in the first 
subsection 4.1.1. In the following subsections, we will go 
beyond scaling arguments and apply functional renormal- 
ization to the FL interaction potential. This will enable us 
to make quantitative predictions for the phase diagram. 

Finally, w e wi ll briefly consider the case with disor- 
der in section 4.2, restricting ourselves again to a scaling 
approach. 



4.1 Pure case Th 



0. 



In the absence of the vdW interaction, the free energy 
density F{£, T) can be represented by the superposition of 
the bare interaction energy z£oKq{1/ X) /P and an entropic 
contribution ^ {T/£'^)/£'^. Here, z is the lattice coordina- 
tion number [z ~ Q for a triangular lattice). The entropic 
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term describes the reduction of the vortex line entropy 
due to its confinement within the cage set up by its neigh- 
bors 1). It can be viewed as a renormalization of the bare 
interaction, taking into account (thermal) fluctuations on 
length scales between A and I. While this renormalization 
of a short range, purely repulsive function can be handled 
within a one parameter renormalization group calculation 

, it is not simply possible to account for the attractive 
tail of the van der Waals interaction within this approach. 

In recent studies, different methods have been applied 
to account for the van der Waals interaction in the Gibbs 
free energy density. In the original work |^ , the bare vdW 
energy, evaluated at the mean FL distance was added 
to the free energy. This means that only contributions 
from the vdW energy on the length scale I are taken into 
account, which leads to a gross underestimation of its in- 
fluence because of its rapid decay for i? > A. In [|oj, on 
the other hand, only the much more important contribu- 
tion on the scale i?min bas been taken into account, where 
-Rmin is the position of the global minimum of the bare 
potential (|8^), see below. This method, which makes use 
of simple scaling arguments, will be reviewed in the fol- 
lowing section. An explicit integration of the fluctuations 
can be performed using a functional renormalization group 
(RG) calculation. The description of this procedure and 
its results is the main focus of this section. The results 
from the RG will be used to check whether the scaling 
arguments lead to a valid description of the effective free 
energy, and to quantify numerical factors that cannot be 
specified within this approach. 

The list of alternative methods is completed by an ap- 
proach that makes use of the mapping of the vortex line 
problem onto 2d-Bosons and a variational procedure. 
This latter approach will be considered in a forthcoming 
publication pa. 



4.1.1 Scaling approach 

The superposition of the bare repulsion and the attractive 
vdW interaction between flux lines, 



V^{R) = 2eoi^o(i?/A) + Kdw(i?), 



(80) 



results in a minimum of their interaction energy at a dis- 
tance -Rniin 

w a\ (<C t) (the typical form of the bare 
interaction potential is shown in the inset of Fig. ||) . 

The position of the minimum can be conveniently quan- 
tified by considering the function w„ (x) = ifo(a:;)— flvdw/s;", 
where n=4 or 5. In Fig. |[ v^(x) is shown for some val- 
ues of flvdw For a huge amplitude range 10^^° < Cvdw < 
1, we find from a numerical fit that the minimum posi- 
tion a„ of Vnix) is well described by the functional form 
a„ ~ Xn + 1.95 ln°'*^(l/avdw) with x/^ = 8.5 and 5:5=12.2. 
From the prefactors of the thermal vdW interaction as 
given in ( [75| ) and ([76|), using the parameters for BiSCCO, 
we read off that the dimensionless vdW amplitude Ovdw 
takes values of order 2 • 10"^ T/K for n = 4 and 10""* T/K 
for 71 = 5. Using the latter value, this corresponds to the 
minimum position ranging from a k, 19 at T = 100 K to 
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Fig. 3. For some values of the vdW amplitude Ovdw, the 
(dimensionless) bare interaction potential vn^ix) = Ki^(x) — 
flvdw/a;* is shown, with a smooth cutoff of the algebraic part 
between x = \ and x = 'h. The exact form of the cutoff will 
be specified below, see Eqs. ( jo^ ) and (^). The inset shows a 
magnified view for two smaller amplitudes flvdw 



a ~ 26 at T = 1 K. Also, the width of the minimmn only 
weakly depends on temperature, being of the order of /5A 
with /3 K, 10. Like a, /? becomes larger for smaller T. 

Due to the strong distance dependence of the vdW 
attraction, its main contribution comes from those con- 
figurations where the line pair is at a distance i?min- In 
the estimate of this contribution developed below, we will 
account for both the thermal and the disorder dominated 
case, which are characterized by the length scales Lt and 
Ldis, respectively. 

We estimate the average vdW interaction by consid- 
ering the configurations of a single line in the absence of 
any FL interaction. With u ~ A(L/LT,dis)^ for the dis- 
placement of a single F L (w ith the roughness exponent ^ 
as discussed in section 2.2), we find from u k, ^ iox the 



mean distance L\\ between two line segments reaching a 
minimum L\\ w LT.dis(^/A)"'^/''. The length Lg of the seg- 
ment over which the line stays in the minimum follows 
from the same argument as ~ iT,dis 0^^^^ ■ Thus, the 
contribution from the vdW attraction to the Gibbs free 



energy density (|79|) is of the order 1 10 
L 



V^dw (^min) 



L 



V^vdw (-Rmin) 



2+i/C 01K 



A2 



(81) 



which is much larger than the vdW interaction at the 
mean distance I. For thermal fluctuations (C — 1/2) the 
mean vdW attraction has therefore the same t dependence 
as the entropic repulsion. In this case, the result can also 
be obtained by mapping the problem onto 2d- Bosons . 

The Gibbs free energy density G[x\ H, T, Tdis = 0) can 
hence be written in the following form ||l^ 



G{x;H, T,0) 



zKo{x) + 



7t — 5t 



— hlnn , 



(82) 
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which has to be minimized with respect to a; = IjX. Ex- 
pression (|8^ ) has to be considered as an interpolation be- 
tween the regimes dominated by the bare interaction at 
high B and the different fluctuation induced interactions 
at low B, respectively. The strength of the entropic repul- 
sion is given by 0] 



7T w 9.08 (T/£oA)' 



(83) 



and the prefactor of the term following from the thermally 
induced vdW interaction is 

0^ T 
a* Eod ln^(7rA/d 

££oA In (Tr/ea) 

The coefficients ct or cannot be determined within 
this scaling approach. This requires taking into account 
contributions from the vdW interaction from all distances, 
according to their statistical weight, instead of singling 
out a the contributions from a typical distance i?min- We 
will quantify the coefficient ct within the functional RG 
calculation presented in the next section. 



4.1.2 Functional Renormalization 

The problem that could not be handled by the approaches 
discussed above is to properly account for the contribu- 
tions from the vdW interaction to the free energy on all 
length scales between A and i. This issue will be addressed 
now by a functional renormalization group calculation. In 
particular, we will use a procedure that is an extension 
of Wilson's approximate recursion relation [^, and that 
has been well established in the context of the wetting 
transition p^ . 

We will only give a short sketch of the procedure here; a 
comprehensive review can be found, e.g., in |29]. The start- 
ing point is the reduced Hamiltonian 7i{r} = H{r}/T 
defined by 



n{r} 



dz 



£o_ 
8T 



dr{z) 
dz 



T 



(85) 



which describes a single line in 1 + d' dimensions, inter- 
acting with a stiff, straight line centered at r = via the 
potential Vb(r). d' is the number of components of r; in 
our case, d' — 2. An ultraviolet cutoff tt/A, correspond- 
ing to the shortest relevant length scale (parallel to the 
FL), is implied. This problem is equivalent to the prob- 
lem of two fluctuating FLs interacting with each other, 
where r(z) = Si(z) — §2(2:) is their difference coordinate. 
Accordingly, the line stiffness eo/4 is equal to the single 
vortex line stiffness ei{k) ~ eo/2 divided by 2, where we 
have approximated the full, dispersive line stiffness (^) by 
its form in the long wavelength limit Afc ^ 1 as given by 



In an RG step, the fast fluctuations r> with wave vec- 
tors between A/b and A are integrated out in an approx- 
imate manner, yielding a renormalized interaction poten- 
tial y(r<). 6 > 1 is the rescaling factor, and r< con- 
sists only of Fourier modes with wave vectors in the range 
O...A/b. Then, one brings back the Hamiltonian to its 
original form by rescaling z ^ z' = z/b and r(z') = 
r<(z = bz') where C is the roughness exponent, as 
defined in (^^. For purely thermal fluctuations which we 
consider here, C = 1/2. For the interaction potential, this 
implies the rescaling 



y(i)[r(z')] - bV'[b'^r{z')]. 



(86) 



This procedure is iterated until all fluctuations on scales 
smaller than £ are integrated out. We denote the renor- 
malized and rescaled potential after the N^^ iteration step 
with FW(r). 

For d' — 2 dimensions, the full nonlinear recursion 
relation is given by 

T/(--)(r).-.51n/^exp(-^-X(r,r')), 

(87) 

where the kernel is defined as 

K(r, r') = — Iv'^^Hb^r - r') + V'^^\b'^r + r') 
2v L 

The length scale a, defined by a? = (r?,), reads 



8T 



duj 1 _ 4T (6 - 1) _ 



A/b 



2n uj^ 



TTEoA 



= a^{b-l), (89) 



and the energy scale v is given by 



v{b) - (1 



with V = TA/ir. 



(90) 



(0), neglecting the contribution 



£2 « 1. 



a and v have been defined such that they do not depend 
on the rescaling factor b, so it is convenient to choose them 
as natural scales for length and energy, respectively. The 
relation of these scales a and v to the physically relevant 
scales A and £0 = LtT/}? takes the following form: For 
the length scale we find — (4/7r^)A^. Here, we have 
identified the short scale cutoff ti/A from ( |85| ) with Lt 
(cf. Eq. 15), the segment length on which the FL makes 
typical perpendicular excursions of order A, because we 
want to integrate out fluctuations of u on scales > A. The 
resulting relation connecting a and A does not depend on 
T. In the following, we will simply set A = a. For the 
energy scale we find v = T^/£oA^ which does depend on 
temperature. 

Note that at low temperatures, Lt, the smallest length 
scale parallel to the FLs in the bare Hamiltonian (|8|) ~ 
which we treat, by virtue of this choice of A, as an effec- 
tive Hamiltonian on the (perpendicular) scale A -, is quite 
large: for T = 1 K it is as large as 2 cm. The applicability 
of the RG is however restricted to samples whose height 
L is much larger than Lt - or, for a given sample, to tem- 
peratures T such that Lt <C L. For smaller systems - or 
smaller temperatures - fiuctuations of the FLs are smaller 
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than A, so it is a good approximation to treat them as 
straight hnes in the dilute hmit H > X we are interested 
in. With L ^ Lt, we are hence in the disentangled flux 
liquid phase where practically no collisions between FLs 
occur [ || . This means that also the contribution from the 
entropic repulsion is absent, hence 77- = in (|8^). Conse- 
quently, it is a good approximation in this limit to replace 
the free energy density F{t\ T) in the Gibbs free energy 
( [79I ) by the bare FL interaction. 



G{l-H,T) 



1 [z 



Vo(^) — Eo/ilnK 



(Lt > L) 

(91) 

where Vq is defined as in (^0|). In the following, we assume 
that Lt ^ L. When discussing the phase diagram, we will 
come back to the opposite case and use expression ( |9l|) for 
the Gibbs free energy density in the low temperature limit 
where Lt ^ L. 

The renormalization procedure sketched above has been 
carried out numerically. Restricting ourselves first to the 
extremely decoupled limit e ^ 0, we define the bare in- 
teraction V'^°\R) = Vo{R)/v as 



V(°n^) = «o ( lMR/\) - avdw/(i?/A)^ 



(92) 



where vq — 2eo/v — 2 (epX/T)^, and the strength of the 
thermal vdW attraction ( [Tq) is determined by Ovdw = 
T/{2deoln^{TTX/dj). f{x) is a function that smoothly cuts 
off the power law tail at i? w A, which we have defined as 




sm TT 



_x-{xi+X2}/2 



X < Xi 

Xi < X < X2 
X > X2 



(93) 

with xi — 1 and X2 — 5. The choice of the cutoff function, 
as well as the actual values of xi and X2, is to some extent 
arbitrary; xi = 1 is however an obvious choice, and X2 
has to be chosen such that the cutoff is not too sharp 
and, on the other hand, does not influence the form of the 
potential in the vicinity of the minimum for those values 
of Cvdw th at are physically meaningful (cf. the discussion 
in section 4.1.1). We have carefully checked that the RG 
results are stable with respect to a variation of X2 ■ 

The potential (|9^), with = 1, is shown in Fig. 3[ For 
BiSCCO, one finds (from the discussion in section 4.1.1 ) 
a typical value Ovdw « 2 • 10"^ at T = 100 K. With this 
van der Waals amplitude, the bare potential Vo{R) has a 
shallow minimum at i?min ~ 18A; cf. the inset of Fig. ^. 

As a consistency check for the procedure, we have com- 
pared the functional renormalization of a purely repulsive, 
short-ranged potential under the action of the recursion 
relation ( p7| ) with a known RG flow, namely the renormal- 
ization of a delta-like, repulsive interaction of two elastic 
polymers by thermal fluctuations. In this model, the am- 
plitude 7 of the delta-like potential vanishes asymptoti- 
cally like 7^(0 ~ 1/^ Here, I is the continuous renor- 
malization parameter which is connected to our N via 
gi _ ^c^'_ niake contact to this model, we have iterated 
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Fig. 4. Functional renormalization of the potential ( |92| ) with 
vo = 10, van der Waals amplitude Ovdw = 3 and a rescaling 
factor 6 = 1.1. Only every 8"^ step N is shown. In the main 
figure, space and energy have been scaled back to original space 
by plotting b-^V'-^\b'^^R). The bare potential with = 
corresponds to the curve with the deepest minimum; for higher 
N, the position of the minimum is shifted to larger values of R, 
and the minimum becomes more shallow, while the amplitude 
in the origin decreases. In the inset, V^^' (R) as defined by (js^ 
is shown. Here, the rightmost curve is the bare potential. 



the recursion relation for the potential ( P4) with Ovdw = 0, 
so that the bare potential is short ranged and purely repul- 
sive, and checked that the amplitude of the renormalized 
potential V'^^^R) does indeed scale like V'-^^O) ~ 1/N 
(the corresponding data for ^'^^^(i?) are not shown here; 
they look however similar to those in the inset of Fig. U, 
see below). 

Now, we have iteratively applied the recursion relation 
( p7| ) to the potential ( p2[ ) for different values of the van der 
Waals amplitude a^dw and global amplitude vq. One has 
some freedom in choosing the rescaling factor b. Larger 
values for b require less steps N to reach the same total 
rescaling factor I = N( In b and hence lead to less numer- 
ical roundoff errors, at the price of a lower resolution in 
I. 

We illustrate the procedure with data based on a bare 
potential with fixed tig = 10, 6 = 1.1 and varying Ovdw 
The evolution of the potential I^'^^(_R) under the action 
of the recursion relation (87) is illustrated in Figs. || and 
H for two different values for a^dw With these values, the 
bare potential Vo{R) belongs to two different generic cases: 
For flvdw < ^vdw where a*,^^ is a critical value that de- 



pends on vq and, very weakly, on b, V^^^ is asymptoti- 
cally mapped onto a purely repulsive short range poten- 
tial, see Fig. ^. For avdw > a*^^, on the other hand, V'^^ 
is mapped onto a purely attractive potential, as shown in 
Fig. 
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Fig. 5. The same plot as in the inset of Fig. ^, now with 
oivdw = 5.7. The potential V''^'(_R) is asymptotically mapped 
onto a purely attractive potential. Again, the rightmost curve 
corresponds to the bare potential. 



4.1.3 Gibbs free energy density 

Now, let us come back to our starting point: The calcula- 
tion of the effective Gibbs free energy density ( fz^ ) which 
we rewrite as 



(94) 



where ^ — —Sohhi k plays the role of a chemical potential. 
G{i; /Lt, T) has to be minimized as a function of £ for fixed 
fi. The effective potential is given by 



fW(A), 



(95) 



where N is determined from £ — Afe*"^. With this N, the 
renormalized V'-^^ includes thermal fluctuations on scales 
between A and the mean flux line distance £. It is eval- 
uated at the (rescaled) position A, corresponding to £ in 
unrescaled coordinates. 

The effective potential is plotted in Fig. |^ for some val- 
ues of Cvdw For Cvdw > aCdwJ where a*^^ « 5.33 with the 
other parameters vq and b set to the values given above, 
V^^{£) has a global minimum at some finite value £min- For 
£ ^ ^min, it decays ~ —l/£'^. This confirms the prediction 
from (^) that the contribution from the vdW interaction 
to the free energy density scales like —l/£'^ in the thermal 
case. 

For Ovdw < a^dw the other hand, V^^{£) is purely 
repulsive and decays for large £ (cf. the inset of 

Fig. revealing that the entropic contribution ~ T/£'^ to 
the Gibbs free energy dominates over the vdW attraction. 

G{£; fjL, T) is now minimized as a function of £ for fixed 
chemical potential /z. The corresponding data for G as a 
function of £ are shown in Fig. ^ for two typical values 
of flvdw smaller and larger than a*^^, respectively. For 
flvdw < ^vdw minimum position £o is continuously 
shifted to larger values when ^ — > from below. This cor- 
responds to a second order phase transition at /i = 0, with 
B = <Po /£q vanishing continuously for H \ H^_^ . Analysis 
of £o as a function of /i reveals B ~ ji ln(l//L(), hence re- 
producing the result from the renormalization group treat- 
ment in [E4|. 
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Fig. 6. Effective potential potential V°^{£) ( psj ) for rescaling 
factor b — 1.1, potential amplitude vo — 10 and several values 
for the van der Waals amplitude avdw in the vicinity of 
The inset shows a double logarithmic plot of a purely repulsive 
¥""{£) with oivdw = 5.0 < fflvdw For comparison, a (dashed) 
line proportional to l/£'^ has been added. 



For Cvdw > the position of the minimum remains 

at a finite position determined by the form of {£) when 
/i approaches from below. This minimum disappears at 
/.t = /i* > 0; for < n < fi*, G has two minima, one 
at a finite length £ and the other one at infinity. This 
scenario describes a first order phase transition with B 
dropping to zero from a finite value By when H \ H^^ = 

Hi - Anfi*/<Po. 

Now, we compare the results from the functional RG 
calculation with the effective Gibbs fre e ener gy density 
derived by scaling arguments in section [l.l.l| . As stated 
above, the scaling behavior for large £ given by expres- 
sion (^) is confirmed by the RG results both in the case 
where the entropic repulsion dominates {St < 7t, corre- 
sponding to Ovdw < Q-vdw) ^^"^ the opposite case where 
the vdW attraction gives the larger contribution. At the 
temperature T* where 7r* = 5t*, one has avdw — o-vdw 
This notion allows for a straightforward calculation of the 
coefficient ct, which determines the amplitude of the vdW 
contribution to the free energy. Using the definitions for 
Vo and Ovdw as given after Eq. (|9^), together with ( |8^ ) 
and (|8^ ) in the limit i?min < d/e, this leads to 



CT ~ 9 



1 



P'^ vq a. 



(96) 



vdw 



vdw numerically for a large range 
10^. Here and below, we will 



We have determined a 
of values for vq — 
report results from functional RG calculations using a 
rescaling factor b — 1.5. It turns out that the right hand 
side of (p^) is indeed basically constant if one defines the 
size of the minimum width to be half the minimum po- 
sition, P = a/2 (which is quite an obvious choice after 
inspection of the form of the bare potential for different 
vdW amplitudes), and uses a = i?min/A with i?niin denot- 
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Fig. 7. Main figure: Gibbs free energy density G{£; fj., T) ( p4| ) 
as a function of the mean flux line spacing £ for different values 
H, based on functional RG data starting from a bare potential 
with Ovdw = 5.6 > a^^^, vo — 10 and b = 1.1. The 
effective potential V'^^{£) for this vdW amplitude is shown in 
Fig. 1^. Only data from the regime < fi < fi* where G has two 
minima are shown. The top curve corresponds to the largest 
value of fi. Inset: The same plot, now for data from a bare 
potential with Ovdw = 5.0 < flvdw The different curves have 
been individually scaled for the plot. From left to right, the 
absolute value of fi becomes smaller. 



ing the true minimum position for a given value of a*^^ 
(for the given range of Vq, a varies between 12 and 19). 
In this way, we obtain ct = 3.9 ± 0.3. The data for the 
critical amplitudes of the vdW interaction a*^^^, as well 
as the values for ct calculated from these data, are shown 
in Fig. ^ as a function of vq . 

Accordingly, we have determined the prefactor in 
the case of finite anisotropy £ > where i?miii > d/e, 
corresponding to the lower term in (p4[). We find a value 
of c'j, = 1.9 ±0.2. 

Now, we will discuss the consequences for the low field 
phase diagram of layered superconductors. We restrict 
ourselves here to an analysis based on the form (|7^) of 
the thermal vdW interaction which decays with the fourth 
power of R and, consequently, to the upper definition of 
St in (p^. We have also calculated the phase diagram ac- 
cordingly using expression (|7^) for the vdW energy which 
is valid for R > d/e. Since the qualitative features are the 
same and even the numerical values are very close to each 
other in the two different cases, we only demonstrate the 
data for the first case. 

We have used the parameters for BiSCCO as given 
above (which implies Ovdw ~ 2 • 10"^ T/K), and have de- 
termined the position Iq of the minimum of the Gibbs 
free energy G{£; ^ — 0,T) which has been calculated from 
the functional RG as described above. In Fig. ^, the re- 
sulting data for the magnetic induction By = <^o/^o ~ 
500G/(^o/A)^ 0] are shown; the corresponding line is la- 
beled sr'''. 
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Fig. 8. 'Critical' van der Waals amplitude Svdw ^ ^ function 
of Vo — 2(eoA/r)^. The product voa*^^ is not a constant; 
instead, a*^^ ^ ^o^'^^ range shown in this figure. The 

values for ct calculated from these data are shown in the inset. 
The straight line corresponds to ct ~ 3.9. 



For comparison, we have determined the same field By 
by minimizing the Gibbs free energy ( |82|) as derived from 
the scaling argument, using (83) and ( ]84| ) for the param- 
eters 7t and St, with ct — 3.9. These oata are shown in 
the same figure, denoted as By'^^ . While the resulting val- 
ues are of the same order of magnitude as those obtained 
from the RG, Bi"'^^ , there is considerable quantitative dis- 
agreement. Hence, although the scaling approach to the 
Gibbs free energy density G{£; /i, T) leads to the correct 
asymptotic scaling of the various contributions, and even 
yields estimates for By in the correct order of magnitude, 
it seems to fail for the exact determination of the mini- 
mum of the Gibbs free energy. Evidently, also the func- 
tional RG approach is not exact: The derivation of ( |87| ) 
contains some severe approximations, and different choices 
of the rescaling factor b lead to slightly different numerical 
results. Nevertheless, the RG represents the most natural 
and adequate method for the calculation of the effective 
Gibbs free energy. 



As discussed above (in section 4.1.2), the RG approach 
makes sense only for temperatures such that Lt is smaller 
than the sample height L. For lower temperatures, corre- 
sponding to the disentangled flux liquid phase, the bare 
interaction enters the Gibbs free energy, cf. Eq. iM). The 
jump in the magnetic induction that results fronitnis ex- 
pression is shown in Fig. ^, labeled Bi^'^^°\ In contrast 
to the RG result (which grows for T 0), 

drops to in this limit. The exact temperature beyond 
which the RG result applies will depend on the sample 
height L, and in a large temperature range around this 
value, the real line By will lie somewhere between B!;^'^* 
andB?"'^^ 

The last line we have included in the phase diagram 
is the lower (first order) melting line Bm{T). This line 
separates the flux liquid phase for B < Bm from a solid 
phase for B > Bm, where the FLs are organized in an 
Abrikosov lattice. For layered superconductors, where the 
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Fig. 9. Low field phase diagram for a layered superconduc- 
tor; physical parameters typical for BiSCCO have been used. 
At low temperatures the Meissner-Ochsenfeld phase becomes 
unstable towards a flux liquid by a first order transition where 
B jumps to Si^'^' , which has been determined from the func- 
tional RG data for the effective interaction energy. For com- 
parison, the same line Bi^'^-' as determined from the Gibbs free 
energy ( |82[ ) that was derived from scaling arguments are shown 
as a dotted line. For very low temperatures, instead of Si'^'^' 
the (dashed) line Bi'"^^'^^ will correctly describe the jump in the 
magnetic induction. Finally, Bm [dashed- dotted line) denotes 
a first order melting transition between the liquid and a solid 
phase at higher magnetic fields B. 



electromagnetic interactions dominates over the Joseph- 
son coupUng, this hne is given by 



BrniT) 



^0 1 -2 



47r cl 
(37r)i/4 



£oA 



(97) 



where cl « 0.3 is the Lindemann number ||30|| . Note that 
Bm > By^'^-' for a large temperature range; only for very 
low temperatures, the two lines cross. For samples with 
finite height L, however, for these low temperatures B^ 
will cross over to the line Bi^'^^'^^ < Bm, so that it is likely 
that By < Bm in the whole temperature range. 

Finally, we want to compare our phase diagram with 
that proposed by Blatter and Geshkenbein [p[|25[|. While 
both proposals basically agree, including the finding that 
By is of order 1 G for low temperatures, there are some 
differences. In particular, these authors find a bubble-like 
shape of the instability region which is bounded by two 
lines By and Bg, where both By and Be are finite for some 
temperature range. In terms of the Gibbs free energy, this 
corresponds to the simultaneous existence of two minima 
at finite values £, which we don't find in the effective po- 
tential as derived from the functional RG procedure, as 
can be seen for example in Fig. 



4.2 Phase diagram in the presence of disorder 

While the functional RG presented in the preceding sec- 
tions only works for the thermal case, the scaling approach 



to th e Gib bs free energy density G{£; H, T, Tdis) from sec- 
tion [4.1.l| can easily be extended to the case with dis- 
order. The contribution from the disorder induced vdW 
attraction can be estimated in the same way as that from 
the thermally induced vdW attraction, which leads to 
Eq. (81), where the value ( = 5/8 for the roughness expo- 
nent has to be used instead of the thermal value 1/2. In 
contrast to the thermal case, where the entropic repulsion 
and the effective vdW attraction have the same depen- 
dence on the mean distance £, the resulting contribution 
from the disorder induced vdW attraction decays faster 
than the (disorder induced) steric repulsion. 

The resulting expression for the Gibbs free energy den- 
sity reads [fol 



G{x;H,T,Tdi,)^G{x;H,T,0) 



7dii 



3/5 



(98) 

The amplitudes of the contribution from the steric repul- 
sion and from the vdW attraction are given by 



7dis ~ Cdis k'*/^ (Tdis/eoA)' 



(99) 



and, for i?min S> d/e, 



Cdis 



3v 



ri-l 



(100) 

respectively, where a and (3 are again the minimum po- 
sition and width of the bare potential, divided by A, and 
the exponent rj was defined in (p9|). 

We briefly recall the essential features of the result- 
ing phase diagram which has been considered semi-quan- 
titatively in Ref. ||l^. The phase diagram in the presence 
of impurities is expected to differ significantly in the low 
temperature region T < Tdis where the disorder mediated 
fluctuations dominate the thermal fluctuations. At T = 0, 
as the magnetic field is increased for weak disorder (but 
A > 1) there is a continuous transition from the Meiss- 
ner phase to the low density phase followed by a jump 
to a high density phase of FLs as the magnetic field is 
further increased. This behavior smoothly crosses over to 
the known thermal phase diagram. As the strength of the 
disorder increases the disorder induced repulsive interac- 
tion dominates and the first order transition at T = 
slowly disappears but still persists at finite low tempera- 
tures. For even stronger disorder the first order transition 
completely vanishes. 

In order to determine the dimensionless coefficients Cdis 
and Cdis, in principal a renormalization procedure simi- 
lar to that applied in the thermal case can be used. Due 
to the presence of disorder and, hence, the existence of 
metastable states, however, this task is considerably more 
complicated. In a closely related model, namely the disorder- 
induced unbinding of a fiux line from an extended defect, 
a functional RG calculation has been applied in it 
would be interesting to see whether the method that was 
developed there can be successfully applied in the present 
case. 
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A Probability density puiSe) 

Consider two independent sets Si, each consisting of N 
Gaussian random numbers £m , with i=l,2. The random 
numbers obey the Gaussian distribution 



B List of symbols 



fie) 



1 



2TTe 



(101) 



Furthermore, define F{e) = de' f{s') which, in this 
case of a Gaussian distribution, is an error function. For 
a given realization, let be the maximal number in Si. 
(Note that by symmetry, the following derivation also ap- 
plies to the minima.) The probability distribution for en 
is given by [53 



gN{eN)^NF'^-\eN)f{eN). 



(102) 



Let pi2{Se) be the probability that the difference between 
the maximal values in 5*1 and S2 is equal to e^' — 
It is given by the convolution 



'N 



Pi2[5e) 



de gN{e)gN{£ + 



(103) 



This integral can be evaluated by a saddle point approxi- 
mation, which results in 



Pi2(<5e) 



In^ 



exp 



In N 5e' 



(104) 



This Gaussian approximation is good for values in the 
range \5e\ < e, which is the regime we are mainly in- 
terested in. For de ^ e, on the other hand, pi2{5e) ~ 
N f(^N + 5e), where ejq = J deegN{e). For the Gaussian 
(pM|), eW^ce In^/^ N with c w 20F/e. 



^vdw 

B 

By 

Ct 

d 

St 

Ao 
A{x) 

A 

e 
£0 
£1 



epin(u) 

e 

It 

G 

H 

K 

£ 
A 
L 
Lt 

Ldis 
^0 



U12 



Vo 

\ 

(...) 



reduced vdW amplitude 
magnetic induction B — 'Po/i'^ 
lattice melting field 
van der Waals transition field 
dimensionless factor in the vdW 
contribution to the free energy 
connected correlation function 
disconnected correlation function 
layer distance 

amplitude of the vdW contribution 
to the free energy 
reduced disorder strength 
wave vector dependent reduced dis- 
order strength 

reduced disorder strength for pan- 
cake vortices A — A{7TX/d) 
anisotropy parameter 
basic energy scale Eq — (^o/47rA)^ 
line stiffness of a single FL 
random pinning potential 
pancake vortex pinning potential 
disorder scale 

amplitude of entropic repulsion 
Gibbs free energy density 
external magnetic field 
lower critical field 
GL parameter k = X/^ 
mean distance between FLs 
London penetration depth 
total flux line length 
thermal length scale Lt ~ e 
disorder length scale 
'chemical potential' fi = —eohhiK 
unit flux quantum — hc/2e 
disorder temperature 
vortex displacement 
typical pancake vortex fluctuations 
dipole-dipole interaction 
basic energy scale in the RG 
= 2£o/v, reduced interaction energy 
London potential 
total bare FL interaction 
thermal vdW interaction 
disorder vdW interaction 
coherence length 
thermal average 
disorder average 



dH) 



(84) 
(16) 
(17) 



(82) 

(H 



(55) 
(83) 
(79) 
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(12) 
(94) 
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